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Relationship among Wigner crystal, charge order and Mott insulator is studied by the path- 
integral renormalization group method for two-dimensional lattices with long-range Coulomb in- 
teraction. In contrast to Hartree-Fock results, the solid stability drastically increases with lattice 
commensurability. The transition to liquid occurs at the electron gas parameter r s ~ 2 for the 
filling 7i = 1/2 showing large reduction from r s ~ 35 in the continuum limit. Correct account of 
quantum fluctuations are crucial to understand charge-order stability generally observed only at 
simple fractional fillings and nature of quantum liquids away from them. 



Coulomb interaction between electrons drives various 
types of electron crystallization ranging from Wigner 
lattice, charge order including stripes to Mott insula- 
tor. The Wigner transition in the continuum space was 
studied on two-dimensional (2D) electron systems by 
quantum Monte Carlo calculations and the transition 
point was estimated to be r s ~ 35 at zero tempera- 
ture ||, where r s = r Q /as with the Wigner Seitz ra- 
dius r = l/^/irn e , the Bohr radius as = Aneti 2 /m*e 2 
and the electron density n e . This 2D Wigner transi- 
tion in quantum region was experimentally realized in 
GaAs heterostructure with extremely high mobility and 
the transition at r s ~ 35 was in good agreement with the 
theoretical prediction fS. 

The most dramatic effect of atomic periodic potential 
in crystal not contained in the electron gas is apparently 
the band formation in electronic spectra due to Bloch 
theorem. In this letter, we stress another crucial concept, 
lattice commensurability, generates different dramatic ef- 
fects when the Coulomb interaction is present. In cor- 
related electron systems in solids as in transition metal 
compounds ^, ||] and organic systems || , charge ordering 
phenomena including stripe-type are common when the 
electron filling n satisfies simple fractional number such 
as 1/2 and 1/3. These phenomena have attracted much 
interest in relation to the mechanisms of high-Tc super- 
conductivity in the cuprates and colossal magnetoresis- 
tance in the manganites. However, in those systems, as 
we discuss later, the effective value of r s is usually esti- 
mated to be rather small in the range r s < 10, though 
we have some uncertainty in the evaluation of the effec- 
tive mass m* and the dielectric constant e. At such low 
values of r s , we normally expect the quantum melting 
of electrons in terms of the above electron gas picture. 
When the electron filling deviates from a simple frac- 
tional number, such charge orders indeed usually melts 
to metals or the charge periodicity is pinned at the sim- 
ple fillings with added carriers localized by disorder. It 
implies that a mechanism of stabilizing the solid only at 
simple commensurate filling works very efficiently. 

Another known result is that the Mott insulator at 
integer fillings can be stabilized even at lower effective 
r s . The Mott insulator can be viewed as the extreme 
limit of strong commensurability. A common feature 
of the charge order and the Mott insulator is that they 



are charge imcompressible state resulted from the short- 
ranged part of the Coulomb interaction. 

In this letter, we show that the mechanism for such 
common stabilization of the solid phase only at simple 
fractional fillings is not clear in the Hartree-Fock(HF) 
calculations, whereas the stabilization and melting can 
indeed be understood only when the quantum fluctuation 
effects arc seriously considered. 

To understand the basic mechanism, we employ a 
Hamiltonian given as 

H = H + Hi, 

H a = - ^ Uj(c\ a Cj a + h.c), 

<i,j>,a 

Hi = U ^2 n^nn + - V l3 n. l n ] , (1) 

i i=£j 

where the notations follow the standard one in the Hub- 
bard type models. We take the long-ranged Coulomb 
term as Vij = V/\r%— Tj\, where we ignore possible screen- 
ing effects arising from electrons in other bands. The 
jelium model is employed by assuming a uniform posi- 
tive charge as a background while the effect of ionic pe- 
riodic potential is represented by the lattice with single- 
band electrons located near the Fermi level. The onsite 
interaction U measured from V may depend on the de- 
tailed structure of atomic orbitals. Here we take U/t = 4 
throughout this paper for simplicity. This ratio may be 
more or less the lower bound for the real situations. In 
any case our results in this letter for n < 1/2 do not sen- 
sitively depend on this ratio. We study spinless fermion 
models in addition to spin-1/2 electrons. For spinless 
fcrmions, the U term does not exist. This corresponds to 
the fully polarized ferromagnetic case. We impose the pe- 
riodic boundary condition in the numerical studies. The 
Coulomb interactions between images across the periodic 
boundaries are taken by the Ewald summation. 

Since our purpose in this letter is to understand a 
generic and universal feature of the interplay between 
quantum fluctuations and commensurability, we restrict 
our study to the square lattice with the anisotropy of 
lattice constants a x and Qjy QbS well as t x and t y taken as 
simple as possible to avoid a possible complexity which 
may intervene the clarification of the essence. The trans- 
fer is limited to the nearest-neighbor pairs in x and y 
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FIG. 2: Hartree-Fock result for dH/dr s vs. r s for 32 spinless 
electrons on 32 by 32 lattices. 

directions denoted by t x and t y as in Fig.l. We choose 
the anisotropy d — a y /a x to make the charge ordering 
on a right triangular lattice possible. By this choice, the 
charge order stability is optimized and can be estimated 
on the same ground as the continuum limit. To see a sys- 
tematic dependence on r s , we control t x and t v by a sin- 
gle parameter of the effective mass m* as t z — fi 2 /2m*a^ 
with z — x or y, which correctly reproduces the disper- 
sion in the continuum limit. 

The above Hamiltonian may be rewritten as Hr = 
H Q + Hj,H = -^E<ij>., ff (4^ + h.c.) - 

^jE<iJ>„,a(4cia + h.C.) + ^(d + i)N e ,Hj = 

u Y,i n ii n ii + \[^i7- e | r"-n i in the ener sy unit of 

Rydbergs, where N e is the total electron-number opera- 
tor. In the continuum limit, by the limit of vanishing a x 
and a y and n — > 0, with r s fixed, the above Hamiltonian 
is reduced to that of conventional electron gas. 

We first show results of the standard HF calculation for 
spinless fermions. For the fillings n < 1/2, the first-order 
transitions are clearly visible from the jump of (dH / dr s ) 
e.g. in Fig. | at r sc = 3.45 ±0.05 for n = 1/32 on 32 x 32 
lattice. Note that the jump of (dH/dr s ) indicates the 
level crossing in the ground state as a function of the 
control parameter r s and gives an evidence for the first- 
order transition. The right ordinate is the scale for the 
second Hamiltonian Hr defined by S = d(r^Hn)/dr s — 
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FIG. 4: The Hartree-Fock results for liquid-solid transition 
points of spinless electrons as a function of inverse of the par- 
ticle number N e in the limit n —* 0. 



(d/Trn)dH /dr s . From Fig. [| the electron crystallization 
is indeed identified in the two-body correlation function 
defined by C(ri) = ± J2 rj ( n rj + r^Tj)- 

For fixed electron number N e , the transition point 
is extrapolated to the continuum limit n — > 0. These 
results are further extrapolated to the thermodynamic 
limit N e — * oo as in Fig. |^. This indicates that the sys- 
tem size dependence is small and the result at N e = 8 
is already a good estimate of the thermodynamic limit. 
At n — 1/2, the transition is of continuous type with the 
absence of a jump in (dH/dr s ) while the peak value of 
the Fourier transform of two-body charge correlation di- 
verges as a Bragg peak with increasing system size for 
> r sc . Then the transition point is estimated as 
r sc = 0.9 ±0.1 after finite-size scaling. 
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To understand effects of spin degrees of freedom, spin- 
1/2 electrons are also analyzed at n = 1/4 and 1/2. The 
transitions appear to be of first order at r sc — 2.05 ±0.05 
and continuous at r sc — 1.35 ± 0.15, respectively. The 
orders of transitions are the same as the spinless case 
and the difference of r sc between these two cases shrinks 
with decreasing filling. We discuss this tendency later. 
We conclude that the unrestricted HF approximation 
predicts the crystallization of 2D electrons ranging from 
r sc ~ 1 at n = 1/2 to r sc ~ 3.69 in the continuum limit. 
We note that the filling dependence is rather small. 

To understand how quantum and many-body fluctu- 
ations modifies the HF results, we have calculated the 
ground state of this system by applying recently devel- 
oped PIRG method (j. This algorithm allows us to 
start from and improve the above unrestricted HF so- 
lutions and reach the correct ground state by taking ac- 
count quantum fluctuations in a controlled fashion. By 
following the path integral formalism, the convergence 
to the ground state is obtained after repeated opera- 
tions of exp[— tH] to Slater determinants with small r 
and by the break-up of the kinetic and Coulomb en- 
ergy terms. The long-ranged Coulomb term is rewritten 
by using the Stratonovich-Hubbard transformation which 
allows this term to operate to the Slater determinants, 
though it increases the number of Slater determinants in 
the linear combination of the wave function With 
increasing the dimension of truncated Hilbert space in 
a nonorthogonal basis optimized by this path integral 
operation, an extrapolation method using the variance 
AE = ((H 2 ) - (H) 2 ) I (H) 2 , is employed to reach the true 
ground state. For the extrapolation, we took the num- 
ber of nonorthogonal Slater determinants up to 256. We 
show results when the extrapolation linearly converges 
well to those of the full Hilbert space. 

We took lattice sizes up to 144 sites (12 x 12 lattice) 
with electrons up to 72 at quarter filling n — 1/2 for both 
spinless and spin-1/2 electron systems (Then for 144 sites 
of spin-1/2 system, we took 36 up and down spin elec- 
trons each). For fillings 1/8 and 1/18, we show only on 
the spinless system for 64 sites and 144 sites, respectively 
with 8 electrons. We restricted to the spinless case for 
these low fillings since the difference of r sc from the sys- 
tem with spins is estimated to be small. This is inferred 
from the comparisons of the PIRG with the HF results at 
n = 1/2 together with the HF results at lower filing with 
spins. In fact, the spin effects appear mainly through the 
exchange process. The exchange is scaled by t 2 where t 
is the effective transfer between neighboring electrons. 
At lowering fillings, it forms lower-energy hierarchy as 
compared with that of the charge order, V. Therefore, 
though spin is important in determining the magnetic 
order, this small energy scale hardly changes the solid- 
liquid transition point for n < 1/2. From comparisons for 
different sizes, the system-size dependence of the transi- 
tion point appears to be small at the low fillings because 
the transition becomes first order. Therefore, though our 
systems at the low fillings are rather small, the results are 
expected to be close to the thermodynamic results. 
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FIG. 5: PIRG result for dH/dr s vs. r 3 for 8 spinless electrons 
on 12 by 12 lattices. 




FIG. 6: PIRG result for two-body correlation function for 8 
spinless electrons on 12 by 12 lattice. 



We first show spinless cases. Fig. |5| shows a typical 
example showing the first-order transition in the PIRG 
calculation at n = 1/18 for 8 fermions on 12 by 12 lattice. 
The behavior of jump is similar to the HF results with 
a clear change of the two-body correlation function as 
in Fig. |(| However, this first-order crystallization takes 
place at substantially higher r s than the Hartre-Fock re- 
sult. The transition at n = 1/8 is also first order while 
that is of continuous type again at n — 1/2. The transi- 
tion point is estimated to be r sc = 1.75 ± 0.25, 13.5 ± 0.5 
and 24.5 ± 0.5 for n = 1/2, 1/8 and 1/18, respectively. 
We have also studied spin-1/2 electrons at n — 1/2 and 
obtained a continuous transition at r sc — 2.0± 1.0, which 
is comparable to the spinless case. 
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FIG. 7: Solid-liquid phase boundary in plane of the fillingn 
and r s . The square indicates the Wigner transition point frj. 
The PIRG results (filled and open circles for spinless and spin- 
1/2 fermions, respectively) are compared with the Hartree- 
Fock results (crosses for spinless and triangles for spin- 1/2 
electrons). The solid curves are guides for the eye. The dotted 
curves show contour lines for the ratio of e/m* to the bare 
value being at 1,2,3,4 and 5 from up to down. The inset 
enlarges small r s region. 

In Fig. 0, the phase boundary of the solid and liquid 
phases, r sc , are shown and the HF and PIRG results are 
compared for n = l/l with integer I. At simple fractional 
filling such as n = 1/2, r sc by the HF approximation is 
relatively good. However, it deviates from the PIRG re- 
sults rapidly with increasing I. The ratio of these two 
estimates increases from ~ 1.5 — 2 at n = 1/2 to ~ 10 
in the continuum limit. In terms of the electron density, 
this difference means from 3-4 to 100, since the density is 
scaled by r 2 s . From the PIRG results, we can understand 
why charge orders are commonly observed in organic and 
transition metal compounds at simple fractional fillings 
as 1/2 and 1/3 while it is pinned or melts away from 
such fillings because r sc increases dramatically with in- 
creasing denominator of the irreducible fraction while the 
parameters of many compounds may lie in this range of 
variation. For example, when we assume a layered per- 
ovskite structure with the lattice constant a = AA and 
the effective ratio e/m* being the twice of the bare value, 
n = 1/2, 1/3 and 2/3 are located in the solid region while 
n = 1/4,3/4 and fractions with higher denominators are 



located in the liquid region as we see in Fig. [7| We note 
that the filling control at fixed lattice constants follows a 
contour line of e/m* in Fig. [?]. As we also see in Fig. [j], the 
HF approximation fails in accounting this general exper- 
imental trend. The general trend can be understood only 
by taking account of quantum fluctuations sufficiently. 

Even metallic states near the charge order may inher- 
ently contain such charge-order and stripe proximities 
e.g. formation of quantum and dynamical defect struc- 
ture, in the region where the static solid could be stabi- 
lized on the mean-field level. Competitions with hierar- 
chical structure formation by various levels of fractional 
numbers of fillings become intrinsic properties leading to 
self-organized structure. This may be particularly con- 
spicuous near the simple commensurate order such as 
the Mott insulator, because the mass renormalization in- 
creases effective r s with more continuous character of the 
transition. Reanalyses of mechanisms of high-Tc super- 
conductivity in the cuprates and the colossal magnetore- 
sistance in the manganites from the present viewpoint 
would be desired in future, since they appear in the prox- 
imity of insulator stabilized by the commensurability. 

In summary, we have studied the liquid-solid transition 
of two-dimensional electrons with long-range Coulomb 
interaction. We have shown how the Mott insulator at 
n = 1, charge orders at < n < 1 and the Wigner lat- 
tice at n — > are connected. The charge order transition 
does not sensitively depend on the spin degrees of free- 
dom for n < 1/2. Although the HF transition points 
are rather insensitive to the electron filling and lattice 
commensurability, the PIRG result with quantum fluc- 
tuations taken in a correct way shows rapid increase of 
r s at the transition with decreasing filling. The present 
PIRG results explain why charge orders or stripes are 
commonly observed only at simple fractional fillings such 
as 1/2 and 1/3 but usually are not away from them. We 
stress that such interplay between quantum fluctuations 
and the commensurability may play crucial role also in 
metals (quantum liquids). It would be interesting if the 
interplay can be experimentally studied by microfabri- 
cation of periodic potential to high-mobility 2D electron 
systems. 
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